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Instantons create a non-local interaction between the quarks, which at finite baryon 
density leads to the formation of a scalar diquark condensate and color superconductiv- 
ity. A mean field approach leads to a self-consistent description of the (qq) and (qq) 
condensates and shows the inevitability of a BCS type instability at the Fermi surface. 
The role of the rearrangement of the instanton ensemble for the QCD phase transitions 
is also discussed. 

1. Introduction 

In the last two decades, instantons were shown to give an explanation of wide range of 
non-perturbative phenomena - from the binding of the lightest hadrons and the formation 
of the chiral condensate to the restoration of the chiral symmetry at temperatures of 
about 150 MeV and the formation of quark gluon-plasma |L|]. Although the detailed 
understanding of the QCD vacuum and of phenomena like the color confinement still 
remains elusive, lattice studies tend to show that the instantons are the principal non- 
perturbative configurations in QCD and their role the phenomena mentioned above can 
be directly observed on the lattice ||. Unfortunately, one area where the lattice studies 
are very difficult due to the complex action is the finite baryon density QCD. It seems that 
the instanton- inspired models are some of the very few methods that work in this area, 
and recent studies |3|-|5| have shown an exciting range of new QCD phases and phenomena 
like the color superconductivity. One can try to follow faithfully the instanton approach 
starting directly from the QCD Lagrangian instead of just considering another 4-fermion 
model interaction. Here we use a mean field treatment of the effective interactions for 
the case of three colors N c = 3 and two flavors Nf = 2. We derive the grand canonical 
potential Q at finite baryon chemical potential fi, containing both the chiral (qq) and the 
color diquark (qq) condensates. We use the exact finite [i instanton formfactors and show 
that the gap equations reveal a BCS like instability at the Fermi surface, which leads to a 
diquark condensation. On the other hand there is a critical coupling for the existence or 
the chiral condensate, which leads to the restoration of chiral symmetry at high densities. 

2. Effective interaction and the finite fi formfactors 

The method of deriving an effective 2A r j-fermion interaction starting from the QCD 
action was first developed in || by D. Diakonov and V. Petrov. It consists of several 
steps that are sketched below. All momenta and coordinates are Euclidean and we use 
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the Euclidean 7 matrices. Starting from the QCD partition function, 

Z = J V^V^VAexp(S QCD ) = J VAdet(#))exp(S gauge ) } (1) 

integrate out the fermions to obtain the fermion determinant and separate the low- 
eigenvalue part from the high-eigenvalue part, by introducing an arbitrary scale 

detOp) = Det low Det high . (2) 

Consider the classical gauge configuration to consist of an instanton-anti-instanton en- 
semble. When sufficiently diluted it can be approximated by the sum ansatz 

A= J £A k + 6A. (3) 

fee/,/ 

The only collective effects are produced by the zero mode part of the fermion determinant. 
It consists of 

T IT = J dtd 3 *(f)\(x - z{) p<p r (x - z T ), (4) 

where <ft are the fermion zero modes. They are known also in the case of finite baryon 
chemical potential [|7| 

*• = W" ( s " - 3 f + + (pffl ) ( C08( '"' ) + ; sin( " r) ) e * (5) 

where p is the instanton radius. 

Next, reintroduce the fermions, but this time as free fields. Reproduce the zero mode 
determinant by calculating a Green's function with Nf(N + + AL + N m ) external legs and 
using the Wick's theorem to do the contractions. 

where in the two flavor case 
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9 + = J dzjdttj n [ / d 4 x^ f (x)i yfaix - zj) J d*y4(y - z T )i ffl> f (y)], (7) 

and the integrals are over the collective coordinates of the instantons. To exponentiate 
the fermion terms, apply inverse Laplace Transformation then go to momentum space 
and do the integrations over the center coordinates and the orientation angles of the in- 
stantons and anti-instantons. The formfactor a(u, k; p) in the resulting non-local effective 
Lagrangian is expressed through the Fourier transform of the fermion zero modes. The 
exact formula depends on the modified Bessel Functions Io(pcu ± /2) and K (po; ± /2) and 



their derivatives, where uj^ = yju + (k± p) 2 . The rather lengthy expression shows that 
the formfactor is centered at the Fermi surface and has a singularity there, but remains 
finite. Its real and imaginary parts are shown on Fig. |l|. 

For the mean field treatment below, it is convenient to write the effective interaction 
known as 't Hooft Lagrangian |§ as a sum of the S, T and U channels of the T matrix 



Figure 1. The real and imaginary parts of the formfactor at \i — 0.5/ p. u and k are also 
in units of 1/p. 



in a tree approximation by doing the appropriate Fierz transforms. The meson channels 
(S+T) have color singlet and color octet terms, the scalar ones are attractive, while the 
pseudoscalars are repulsive: 

C me s = [(^r-^) 2 + (^r- 75 ^) 2 ] + 16N f { ~, 2 _ 1} [(^~A» 2 

+ ^T-X^Y] - 32iVe( ^_ 1) (^"^A» 2 | • (8) 

The diquark channels contains both color antisymmetric 3 and symmetric 6 terms, the 
scalar and the tensor are attractive: 



C-diq — G 



1N2(N C -1) 
1 



(iP T CT 2 \ a A ^T 2 \ a A Ci; T ) + (VCr 2 X a A75 i;)^T 2 X A7 ^ T ) 



+ 16N 2 {Nc + 1} ^ T Cr 2 X%a^)^r 2 \ a s a, u C^ , (9) 

where C is the charge conjugation matrix, t 2 is the anti-symmetric Pauli matrix, \a,s 
are the anti-symmetric (color 3) and symmetric (color 6) color generators (normalized in 
an unconventional way, tr(X a X b ) = N c S ab , in order to facilitate the comparison between 
mesons and diquarks). When written in momentum representation, each fermion has a 
formfactor y/a attached to it. The overall coupling constant G comes from the expo- 
nentiation of the fermion terms through the inverse Laplace transformation. One can 
do a saddle point approximation of this integral by varying G and thus relate it to the 
parameters of the instanton ensemble - total density ni nst and average radius p. 
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3. Mean field approach 

There is a variety of mean field methods, that yield equivalent results. One of the 
most convenient is the use of a bilocal effective action T[G,F,F] [PHO]], which depends 
on the propagators G(x,y) = {Tip{x)%[} F(x, y) = (Tip(x)ijj T (y)) , and F(y,x) = 
(Tip T (x)ip(y)). In momentum representation one obtains: 

T[G,F,F] = -l/2trln[G(-p)] - l/2tr ln[G(p) + F(p)(G T (-p))- 1 F(-p)] 

+l/2tr[G 1 (p)G(p) + G 1 (-p)G(-p) - 2] + V[G, F, F]. (10) 

Assuming the breakdown both of the chiral and the color symmetry through scalar color 
singlet ("0V) an d color 3 {ip T C^^T 2 ^ ■ Xa^) condensates, one can write the relevant for the 
mean field approach terms of the potential V in the zeroth meson-diquark loop approxi- 
mation 

V[G, F, F] = Gg^(tr[G! + G 2 ]|a|) 2 + 2 _ (tr[Ag(Gi + G 2 )|«|]) 2 

+Gg^^ ri ytr[C7 5 T 2 A A Fa(p)]tr[C 75 T 2 A A Fa(-p)]. (11) 

The unit vector £ determines the direction of the diquark condensate in the space of 
Xa- Without restrictions we can chose it to point towards A 2 . The breakdown of the 
color symmetry leads to the formation of two different chiral condensates, one involving 
the quarks of the first two colors, and the second - for the quarks of the third color. 
Correspondingly, we split the propagator G = G\ + G 2 , with G\ containing the projector 
P to the upper left 2x2 square of the SU(3) matrices, while G 2 containing the projector to 
the lower right entry. Because both G\ and G 2 mix with the A§ matrix, one has to include 
the color octet interaction term, which contains A§. (Please, note again the non-standard 
normalization of the A matrices.) Varying V one obtains the Hartree-Fock equations for 
the propagators. In a picture representation they are shown in Fig. |2|. 
The solution of these equations is 

-(i j,- + m2 )(l-P) (ifi++ mi )P 

^2 — / xn , n > "1 



(p ) 2 + m 2 (i fa — mi)(i fa + + mi) + A 2 ' 

-iT 2 X 2 C'y 5 Aa(-p) 
(i fa + + mi)(i fa~ — mi) + A 2 



where far* 1 = 7o(cu ± ifj) + jk and mi j2 = m + Mi 2 , with m being the current quark mass 
and Mi >2 = a(u, k, mw)Mi i2 , Ai )2 = a(u), k, mu) Al, 2. Substituting into (|T0l) one obtains 
the grand canonical potential 

n(M u M 2 ,A^) = Aj ^ln[(K + ) 2 + A 2 )(K-) 2 + A 2 )] + 2/ ^ M^) 2 ^ ) 2 ] 
d 4 p m 1 M 1 (s 1 + A 2 - /x 2 ) + A 2 (si + A 2 + /i 2 ) /• rf 4 p m 2 M 2 (s 2 - /r 2 ) 
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(2tt)* (K+) 2 + A 2 )((cr) 2 + A 2 ) •/ (2tt) 4 (c^) 2 ^" 

d 4 p|a| / 2mi(si + A 2 — fi 2 ) m 2 (s 2 — fi 2 ) 



+ 9° 
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Figure 2. Hartree-Fock equations for the propagators. The effective masses Mi,M 2 and 
the gap A are proportional to the corresponding chiral and color condensates. 



d A p\a\ I m 1 (s 1 + A 2 — /i 2 ) 
■(2^0* I (K+) 2 + A 2 )((^D 2 + A 2 ) 



m 2 (s 2 ~ A^ 2 )' 
G4) 2 K) 2 . 



<i 4 p|a| 



A(si + A 2 + /x 2 ) 



(2tt) 4 ((^+) 2 + A 2 )((cur) 2 + A 2 ) 



,(13) 



where (u;f 2 



a; 2 + 



A; 2 + m 2 2 ± /i) 2 , si )2 = cu 2 + A; 2 + m 2 2 . By differentiating with 

respect to Mi, M 2 and A one can obtain a system of coupled gap equations, which can be 
solved numerically for each value of \i. An easier approach for studying the phase diagram 
of the system is to plot Q as a function of two of the parameters for different values of \i 
and a fixed third parameter. The minima will correspond to the different possible phases 
and the phase transitions correspond to two minima having an equal value (first order), 
or joining together (second order). But even without a computer, one can study the 
limiting cases - pure chiral condensate (Mi = M 2 , A = 0), or pure diquark condensate 
(Mi = M 2 = 0). The corresponding gap equations have the following form: 
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[s! + A 2 + fi 2 )\a\ 
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(2vr)MK + ) 2 + A 2 )(K-) 2 + A 2 )' 

d 4 p {S\ — /U 2 )H 2 

(27T) 4 (S1 -fi 2 ) 2 + 4/i 2 W 2 ' 



(14) 
(15) 



The formfactor \a\, which enters in Mi,M 2 , A, although having a cusp at k — \i is also 
continuous, so that for the qualitative picture one can substitute it with a sharp cutoff 
around the Fermi surface. If we remove the gap A the first equation develops a logarithmic 
divergence of the BCS type and can not be satisfied for any G: 



1 

96V 



v+t k 2 dk 

ju-e k + /i 



+ 



v+t k 2 dk 



k 



+ 



k 2 dk 



fj.—e jJi 



(16) 
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The last two terms, which are divergent correspond to quarks and holes close to the Fermi 
surface k = //. The second equation, however, is of Nambu type 

, = * CF, f" +< , k2 " k (17) 

and needs a minimal critical coupling for a solution to exist even at finite \x. Note also 
that the Fermi momentum kp = i//i 2 — m 2 . Similar conclusions were reached in H using 
a different method and ad-hoc formfactors. 



4. Phase Transitions 

Although in the picture drawn so far there can be phase transitions, they may be very 
dependent on the particular scenario, e.g. in || the chiral restoration transition is due 
mainly to the formf actor which is centered not at k — fi, but at k = 0, so that at sufficiently 
large chemical potential, it cuts off effectively the interaction. Another possibility is the 
Debye screening of the instantons in the plasma phase, which also leads to decrease 
in the effective four-fermion coupling. Studies of the finite temperature chiral restoration 
transition [ I^JI3 |, has shown the importance of the correlations in the instanton ensemble, 



which grow with T. One can take into account these correlations in the mean field 
approach by considering the tightly bound and polarized instanton - anti-instanton pairs 
as new objects called molecules, which form a new component. The phase transition 
is characterized by complete pairing of all instantons into molecules. The role of the 
molecules is to deplete the random component, which is responsible for the collectivization 
of the zero modes and it also provides additional four-fermion interactions. 

At finite chemical potential, two mechanisms for suppressing (ipip) may exist: the de- 
pletion of the random component by the molecules, and the competition from (i/ji/j) - i.e. 
instantons previously engaged in the infinite "vacuum quark clusters" (the collectivized 
zero modes) become engaged in the "matter diquark chains" . 

If fi > > m s , we should consider the s quark as a light one and study the phase diagram 
as a function of its mass too. As recent papers has shown Qthe phase structure in the 
Nf = 3 case is richer and more exotic. There is the possibility of many new condensates, 
and patterns of symmetry breaking, i.e. SU(3) C x SU(3)f — > SU(3)d- The color super- 
conductivity is complete, because there is a gap for all quarks. In this case the instantons 
can not form alone a diquark condensate, but the molecules can. Their contribution to 
the diquark binding may be not a correction, but the main one. 



5. Conclusion 

The mean field approach to the instanton induced quark interaction has proved itself 
as a reliable method to study the richness of QCD phenomena at finite baryon density, a 
situation in which very few other methods seem to work. It has shown that the expecta- 
tions of a BCS type of superconductivity transition were correct, and simultaneously they 
show how the chiral symmetry is restored. However, if one wants to make this approach 
more realistic, he should take into account the correlations and the restructuring of the 
instanton ensemble, something that has been shown to play the key role in the finite 
temperature phase transition. 
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